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ABSTRACT 

We present deep observations of the 12 CO (1—0) and (3—2) lines in the ultra-luminous infrared and radio galaxy 4C12.50, 
carried out with the 30 m telescope of the Institut de Radioastronomie Millimetrique. Our observations reveal the cold 
molecular gas component of a warm molecular gas outflow that was previously known from Spitzer space telescope 
data. The 12 CO(3 — 2) profile indicates the presence of absorption at — 950 kms -1 from systemic velocity with a central 
optical depth of 0.22. Its profile is similar to that of the HI absorption that was seen in radio data of this source. A 
potential detection of the 0— >1 absorption enabled us to place an upper limit of 0.03 on its central optical depth, and 
to constrain the excitation temperature of the outflowing CO gas to >65 K assuming that the gas is thermalized. If the 
molecular clouds fully obscure the background millimeter continuum that is emitted by the radio core, the H2 column 
density is >1.8xl0 22 cm -2 . The outflow then carries an estimated cold H2 mass of at least 4.2xl0 3 Mq along the nuclear 
line of sight. This mass will be even higher when integrated over several lines of sight, but if it were to exceed 3xl0 9 Mq, 
the outflow would most likely be seen in emission. Since the ambient cold gas reservoir of 4C12.50 is l.Ox 10 Mq, the 
outflowing-to-ambient mass ratio of the warm gas (37%) could be elevated with respect to that of the cold gas. 
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1. Introduction 

Through jets, winds, and radiation pressure, active galactic 
nuclei (AGN) can affect the collapse of their surrounding 
molecular gas and the formation of new stars in their host 
galaxies. While simulations attribute a considerable role to 
AGN feedback in shaping galaxy properties (e.g., Croton 
et al., 2006; Sijacki et al., 2007; Booth & Schaye, 2009; 
Hopkins & Elvis, 2010; Debuhr et al., 2012), the question 
remains open from an observational point of view. Do the 
AGN-driven outflows that suppress star formation occur 
frequently enough to affect the observed luminosity and 
mass functions of galaxies? To answer this question, we 
not only need to identify objects with massive, AGN-driven 
outflows of molecular gas, but also to take into account the 
outflow effects on the different phases of this gas. 

The various phases of an AGN-driven molecular gas 
outflow have been studied in a nearby ultraluminous in- 
frared galaxy (ULIRG), Mrk231. Plateau de Bure inter- 
ferometric observations showed that the outflowing CO, 
with a corresponding H2 mass of 6xl0 8 M Q and a flow 
rate of 700M Q yr _1 , is capable of suppressing star forma- 
tion (Feruglio et al., 2010). Aalto et al. (2012) discovered 
the outflow's dense (>10 4 cm~ 3 ) component from broad 
HCN(1— 0) profile wings. On the basis of profile fitting of 
the 79 /im and 199 /im OH lines seen with Herschel (Fischer 
et al., 2010), the outflow contains >7x 10 8 M of molecules. 
In total, more than 10 9 M Q of gas can be in the wind. 

Another straightforward and quantitative comparison 
that has yet to be performed is that between the warm 
and the cold H2 gas mass in the outflow vs. in the am- 
bient interstellar medium (ISM). A source that is suitable 
for this analysis is 4C12.50. It is the only source known 



to have a massive outflow of gas of few hundred Kelvin, 
as inferred from its mid-infrared H2 rotational line pro- 
files (Dasyra & Combes, 2011). In this letter, we present 
evidence of cold molecules in its outflow, and we evaluate 
the relative fractions of warm and cold gas that the AGN 
feedback kinematically distorts. We adopt Ho=70 kms -1 
Mpc~\ ft M =0.3, and' ft A =0.7 throughout. 



2. Data acquisition and reduction 

The observations were carried out with the 30 m tele- 
scope of the Institut de Radioastronomie Millimetrique on 
December 21-22 2011, and January 1 2012 (as part of the 
program 235-11). For the 12 CO(l - 0) observations, EMIR 
was tuned to 102.700GHz. For the 12 CO(3 - 2) observa- 
tions, the receivers were tuned to 308.086 GHz during the 
first observing run and to 308.600 GHz during the second 
observing run, to ensure that the observed line properties 
were unaffected by possible standing waves, atmospheric 
line residuals, or bad channels. Parallel to the CO observa- 
tions, we carried out HCN and HCO + (2—1) observations 
with a common tuning of 158.427GHz. Both the FTS and 
WILMA backends, with corresponding resolutions of 0.2 
MHz and 2MHz, were simultaneously used as backends to 
allow consistency checks. The inner and outer parts of the 
FTS were used, enabling us to better sample the continuum 
near each line. The wobbler switching mode with a throw 
of 80" and 0.65 sec per phase was used to rapidly sample 
the continuum. The telescope pointing accuracy was 2", 
and the system temperature varied from 310 K to 580 K at 
1 mm. 
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Fig. 1. Left: FTS and WILMA spectra of 12 CO(3 - 2) in 4C12.50, binned into 62 kms" 1 channels and shifted by a constant Co, 
when needed. Right: Flux-calibrated line profile of 12 CO(3 — 2), averaged over the WILMA and the FTS data to be least affected 
by potential artifacts. Gaussian functions with parameters that most closely fit the emission component at the systemic velocity 
and the absorption component at — 950 kms -1 are shown in blue and cyan, respectively. The inverse HI absorption that is seen 
in the radio data of 4C12.50 (Morganti et al., 2004) is overplotted in orange for an arbitrary scale and continuum level. 



To reduce the 1 mm data, we removed bad channels from 
all individual scans. The three parts of each FTS scan, cor- 
responding to the output of a different spectrometer, were 
fitted with linear baselines and brought to the same con- 
tinuum level, i.e., to the median value of the full scan. This 
step was performed to ensure that there is no continuum 
discontinuity between the three parts that could introduce 
artificial features into the combined spectrum. We averaged 
the individual scans of each frequency tuning and back- 
end, discarding scans with unstable baselines. The inner 
and the outer FTS spectra were then stitched together. For 
each backend, the spectra of the two frequency tunings were 
brought to a common reference frequency of 308.245 GHz 
and averaged. The final spectra are shown in Figure 1. Their 
total on-source integration time is 2.7 hours. 

At 3 mm, we averaged all useful WILMA or FTS scans, 
after shifting them to a reference frequency of 102.756 GHz 
and removing bad channels from them. In contrast to 
the 1 mm data, the 3 mm data had signatures of stand- 
ing waves owing to the strong synchrotron continuum ra- 
diation, which increases with A. We iteratively fitted the 
baseline of the WILMA data with sinusoidal functions of 
periods that are integer multiples of each other, represent- 
ing different harmonics of the same standing wave. The 
best-fit solution that was indicated by the CLASS baseline 
routine corresponded to functions with frequencies of 0.28, 
0.84, 1.68, and 2.52 GHz, and an amplitude of lmK. These 
functions were removed from both the WILMA spectrum 
and the central part of the inner FTS spectrum. The fi- 
nal 12 CO(l — 0) spectra, with a total on-source integration 
time of 1.2 hours, are shown in Figure 2. A similar approach 
was taken for the 2 mm data. A main-beam-temperature- 
to-flux conversion factor of 5.0 JyK -1 was used to flux 
calibrate the final spectra at all wavelengths. 



3. Results: CO in emission and in absorption 

The 12 CO(3 - 2) profile of 4C12.50 consists of an emis- 
sion component at the galaxy systemic velocity, i.e., at 
z=0.1218, and of an absorption feature at approximately 
— 1000 kms -1 from it. The CO(2— >3) absorption is seen in 
all four frequency tuning and backend combinations, indi- 
cating that the feature is intrinsic to the source. Its signal- 
to-noise ratio (S/N) is ~6 in both the WILMA and the FTS 
data, which have a corresponding root-mean-square noise of 
1.7mK and 1.8 mK, respectively. The outer FTS data pro- 
vide further evidence that the absorption is real by ruling 
out the existence of standing waves that could create a fea- 
ture of this depth. The observed optical depth at the center 
of the line, r o b s ,0j is 0.22 (±0.05), with an absorption min- 
imum at 30 (±2) mJy below the local 150(±20)mJy con- 
tinuum. Averaging the best-fit Gaussian parameters of the 
data for each backend indicates that the minimum is at 
—950 (±90) kms -1 . The CO absorption profile remarkably 
resembles that of the HI absorption seen in radio data 
(Figure 1; Morganti et al., 2004). Its resolution-corrected 
width is 250 (±80) kms -1 , reflecting the collective motions 
of multiple clouds along the line of sight. 

An absorption line is also seen at —1100 kms -1 away 
from 12 CO(l — 0) after the progressive subtraction of 
standing wave harmonics from the 3 mm baseline. However, 
the wave amplitude is comparable to the absorption mini- 
mum (Figure 2), impeding its exact profile study. We thus 
simply constrain r b s ,o(0— sT) towards the 0.51 (±0.06) Jy 
continuum background to an upper limit of 0.03. 

The 12 CO(l — 0) emission peaks at a main beam tem- 
perature of 9.2 (±1.2) mK above its underlying contin- 
uum. Its resolution-corrected width is 400 (±95) kms -1 , 
and its intensity, ico(i-o); is 4.0 (±1.3) Kkms -1 (see also 
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Fig. 2. WILMA (left) and FTS (right) 12 CO(l - 0) spectrum of 4C12.50, binned to 93 kms" 1 channels. The dashed histogram 
corresponds to the average spectrum prior to the removal of sinusoidal baselines, which represent the harmonics of a standing wave 
(shown in orange). The final spectrum is shown as a solid, filled histogram. 



Evans et al. 1999, 2005). From this, we deduce that the 
cold H 2 gas reservoir of 4C12.50 is 1.0 (±0.3)xl0 10 M Q , 
assuming that the gas mass is given by the product 
23.5 a I C o(i-o)Vb D 2 L /(l + z) 3 M Q (Solomon et al., 1997), 
where fls is the telescope main beam area in arcsec 2 , Dl is 
the source luminosity distance in Mpc, and a is the CO lu- 
minosity to H2 mass conversion factor. We used a 103 GHz 
beam of 28" (61kpc), and a=0.8M Q (Kkms" 1 pc 2 )" 1 
(Dowries & Solomon, 1998) because 4C12.50 has a 
ULIRG-like infrared luminosity (2.5xlO 12 L0) based on its 
12, 25, 60, and 100 ^m flux (Golombek et al., 1988; Sanders 
& Mirabel, 1996; Guillard et al., 2012). 

The 12 CO(3 - 2) emission is detected with a S/N of 
25 and 21 in the WILMA and the FTS data, respec- 
tively. Its main beam temperature peaks at 17.5 (±1.0) mK 
above the local continuum. Its resolution-corrected width is 
510 (±65) kms -1 , and its intensity is 9.9 (±1.6) Kkms -1 . 

The /co(3-2)/^co(i-o) ratio of the ambient gas, 2.5, is 
below its theoretically predicted value of 9 for the case of 
optically thick gas that is thermalized to a single T cx . It 
is also below that of the central line of sight of M82, of 
several local AGN, ULIRGs, and radio galaxies, but above 
that of the outflow of M82, and similar to that of the Milky 
Way. This comparison is visualized in Figure 3 for the line 
luminosities. Even though this ratio is indicative of large 
reservoirs of cold and diffuse gas, it does not rule out the 
presence of a warm and dense, highly excited CO com- 
ponent, which could be detected through high J emission 
lines. The non-detection of HCN (or HCO + ) at the 3a level 
of 0.86 K kms" 1 for a width of 460 kms" 1 cannot be used 
to argue for or against either interpretation. It leads to 
Ico/Ihcn>8 (Gao & Solomon, 2004) for a linearly interpo- 
lated CO(2-l) intensity of 7 Kkms" 1 . 

4. Discussion: on the properties of a multi-phase 
AGN-driven outflow 

The rapid molecular gas motions that we detected pro- 
vide direct evidence of an outflow in 4C12.50 because none 
of its merging components (Axon et al., 2000) is moving 



faster than ±250 kms" 1 from its systemic velocity (Holt 
et al., 2003; Rodriguez Zaurfn et al., 2007). The maxi- 
mum gas velocity reaches —1500 kms" 1 , which is atypical 
of supernova-driven outflows (e.g., Heckman et al., 2000; 
Rupke et al., 2005), but characteristic of AGN-driven out- 
flows (e.g., Rupke & Veilleux, 2011; Sturm et al., 2011). 
AGN radiation pressure and winds can affect the gas kine- 
matics in 4C12.50: the AGN radiation does reach and ionize 
part of the outflowing gas, which contains both Ne V and 
OIV ions (Spoon & Holt, 2009; Dasyra et al., 2011). The 
gas can also be pushed by shocks created as the radio jet 
of 4C12.50 (Stanghellini et al., 1997) propagated through 
the ISM. An HI outflow in front of a bright radio knot 
(Morganti et al., 2004) makes this scenario plausible for the 
neutral gas, and even for the molecular gas. A jet-ISM in- 
teraction can be responsible for the irregular H2 kinematics 
(Dasyra & Combes, 2011; Nesvadba et al., 2011), the strong 
H2 emission (Ogle et al., 2010), and the high CO excitation 
(Papadopoulos et al., 2008) of radio galaxies. 

To constrain the line-of-sight CO column density, Nco, 
from the observed absorption, we assume that the gas is ho- 
mogeneously distributed and in local thermodynamic equi- 
librium (LTE). We can then relate A*co to r through 
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(e.g., Wiklind & Combes, 1995), where c is the lightspced, 
h is the Planck constant, k is the Boltzmann constant, Ej 
and gj are the energy and the statistical weight of the state 
J, respectively, Aj+\ t j is the Einstein coefficient for spon- 
taneous emission, v is the frequency of the emitted line, 
and Q is the partition function at the excitation temper- 
ature T ox . The ratio / r(2 -> 3) dV / J r(0 ->■ 1) dV agrees 
with its observed limit, >6.4, if T CX >65K. Turbulence and 
shocks can easily bring the gas to this T cx through colli- 
sions. For T CX =65K, A^co is equal to 1.8xl0 18 cm" 2 . 

This number will increase with increasing T cx , but will 
remain on the same order of magnitude up to 200 K. The 
filling factor of the background source by the molecular 
clouds, / c , can play an important role in the T ox and 
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iVco determinations. Partial coverage of the background 
source (/ C <1) can cause the observed (beam-averaged) op- 
tical depth to deviate from its actual value, which is then 
equal to -;n[l-(l-e- r ° b =)// c ] (Wiklind & Combes, 1994). 
For gas in LTE, this relation is computed from the radiative 
transfer equation, assuming that T cx <<Tb g and that the 
cosmic microwave background contribution to the gas ex- 
citation is negligible (e.g., Staguhn et al., 1997). Requiring 
the gas to be optically thick limits f c to values above 0.2. 
Data of multiple transitions are needed to measure / c , fur- 
ther constrain Nco j & n d examine any large velocity gradi- 
ents (Scoville & Solomon, 1974) and different energy level 
populations (Israel et al., 1991). 

For the present assumptions and a CO-to-H2 abundance 
ratio of 10 , the H2 column density towards the out- 
flow, iV H2 , is >2xl0 22 cm" 2 . The HI column density is 
3xl0 21 cm" 2 , based on both optical Nal data (Rupke et 
al., 2005) and radio HI data (Morganti et al., 2005, for 
T sp i„ = 1000K). For the ambient HI gas, which is located 
100 pc northwest of the radio core, the column density is 
6xl0 20 cm" 2 -10 22 cm" 2 (Morganti et al., 2004; Curran & 
Whiting, 2010, for T sp ; n = 100K). The largest uncertainty 
in this comparison is the wavelength-dependent spatial dis- 
tribution of the background and the foreground source. The 
radio core emission typically dominates the millimeter and 
even the submillimeter (850 /im; Clements et al. 2010) con- 
tinuum, because the emission from diffuse jet knots is de- 
scribed by a steep power law that declines rapidly with 
frequency (Krichbaum et al., 1998, 2008). We thus assume 
that the millimeter background source area is contained 
within the 2 cm beam. 

For a beam area of 3 milliarcsec 2 (Lister et al., 2003) 
and A H2 =1.8xl0 22 cm" 2 , the outflow entrains 4.2 xlO 3 M Q 
of cold H2 gas in front of the millimeter core. How this 
mass compares with the mass of the cold H2 gas that flows 
over all lines of sight is unknown. It therefore remains un- 
clear whether the outflow carries less cold than warm H2 
gas. However, the outflowing-to-ambient mass ratio of the 
cold gas will be lower than that of the 400 K gas (37%; 
Dasyra & Combes 2011), as long as the cold outflowing gas 
weighs less than 3xl0 9 Mq. If this amount of gas existed 
in the millimeter beam, and without a background source 
illuminating most of the clouds, the outflow would be seen 
in emission. Our 12 CO(0 — > 1) data instead limit any ad- 
ditional outflowing cold gas mass to be <2xl0 9 M Q when 
integrating over a 1000 kpc 2 area and a 1500kms _1 veloc- 
ity range (Batcheldor et al., 2007). 



5. Summary 

We have presented deep 30m observations of 12 CO (1—0) 
and (3—2) in the nearby ULIRG and radio galaxy 4C12.50, 
which has an impressive outflow of warm (400 K) H2 gas 
(Dasyra & Combes, 2011). Our observations led to the dis- 
covery of the cold component of the molecular outflow. 
Absorption from ,1=2—5-3 is seen in data obtained with 
two different frequency tunings and backends. It peaks at 
— 950kms _1 , and its profile resembles that of the HI ab- 
sorption component that traces the neutral gas outflow 
(Morganti et al., 2004). Absorption is also potentially seen 
for the 0— s-1 transition. Treating the optical depth of the 
latter as an upper limit and assuming that the clouds 
fully cover their background source, we have found that 
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Fig. 3. CO line luminosities in 4C12.50 and other sources, nor- 
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Fixsen et al. (1999), Dumke et al. (2001), Weiss et al. (2001, 
2005), Papadopoulos et al. (2007, 2008), Panuzzo et al. (2010), 
van der Werf et al. (2010), and Krips et al. (2011). Typical un- 
certainties are 0.1-0.3 dex. 



the molecular gas has an excitation temperature of at 
least 65 K in LTE conditions. For this temperature, A^h 2 
is 1.8xl0 22 cm" 2 . This suggests that the outflow carries 
4.2xlO 3 M of cold H2 gas along the nuclear line of sight. 
Even if the mass of the cold outflowing gas could exceed 
that of the warm outflowing gas, 5.2xl0 7 M Q , for a higher 
Nco value or when integrating over several lines of sight, it 
would most likely not exceed 3xlO 9 M , i.e., one third of 
the ambient cold gas reservoir. The outflowing-to-ambient 
mass ratio could thus be elevated into the warm gas phase 
with respect to the cold gas phase. 
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